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In contrast to all other debris disks, where the dust can be seen via an infrared excess over the stellar photosphere, the dust emission of 
the Edgeworth-Kuiper belt (EKB) eludes remote detection due to the strong foreground emission of the zodiacal cloud. In this paper, 
we access the expected EKB dust disk properties by modeling. We treat the debiased population of the known transneptunian objects 
(TNOs) as parent bodies and generate the dust with our collisional code. The resulting dust distributions are modified to take into 
account the influence of gravitational scattering and resonance trapping by planets on migrating dust grains as well as the effect of 
sublimation. A difficulty with the modeling is that the amount and distribution of dust are largely determined by sub-kilometer-sized 
bodies. These are directly unobservable, and their properties cannot be accessed by collisional modeling, because objects larger than 
(10 . . . 60) m in the present-day EKB are not in a collisional equilibrium. To place additional constraints, we use in-situ measurements 
of the New Horizons spacecraft within 20 AU. We show that, to sustain a dust disk consistent with these measurements, the TNO 
population has to have a break in the size distribution at j < 70 km. However, even this still leaves us with several models that all 
correctly reproduce a nearly constant dust impact rates in the region of giant planet orbits and do not violate the constraints from 
the non-detection of the EKB dust thermal emission by the COBE spacecraft. The modeled EKB dust disks, which conform to the 
observational constraints, can either be transport-dominated or intermediate between the transport-dominated and collision-dominated 
regime. The in-plane optical depth of such disks is T||(r > 10 AU) ~ 10"* and their fractional luminosity is fj ~ 10"^. Planets and 
sublimation are found to have little effect on dust impact fluxes and dust thermal emission. The spectral energy distribution of an 
EKB analog, as would be seen from 10 pc distance, peaks at wavelengths of (40 . . . 50) \im at F « 0.5 mjy, which is less than 1% 
of the photospheric flux at those wavelengths. Therefore, EKB analogs cannot be detected with present-day instruments such as 
Herschel/PACS. 
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1. Introduction 

The Edgeworth-Kuiper Belt (EKB) with its presumed collisional 
debris is the main reservoir of small bodies and dust in the 
Solar System and constitutes the most prominent part of the 
Solar System's debris disk. However, the EKB dust has not been 
unambiguously detected so far The observational evidence for 
the EKB dust is limited to scarce in-situ detections of dust in 
the outer Solar System by a few spacecraft, pa r tly with uncal- 
ibrated "chance detectors " dGurnett et al.lll997t iLanderaf et al.l 
l2002t iPoppe etalJbol O^. In addition, there are rough upper 
limits on the amount of dust from the non-detection of ther- 
mal emi ssion of the EKB du st on a bright zodiacal light fore- 
ground (iBackman et al.l[T995 b. Given the lack of observational 
data, one can only access the properties of the EKB dust by 
modeling. Such a modeling takes the known EKB populations 
to be parent bodies for dust and uses collis i onal models to gen- 
erate dust distributions (ISternlll995l Il996t I Vitense et al1l2010t 

lKuchner&Stari3l2010h . 

In our previous paper (IVitense et al.ll2010h . we took the cur- 
rent database of known transneptunian objects (TNOs) and em- 
ployed a new algorithm to eliminate the inclination and the dis- 
tance selection effects in the known TNO populations and de- 
rived expected parameters of the "true" EKB. Treating the debi- 
ased populations of EKB objects as dust parent bodies, we then 
produced their dust disk with our collisional code. 
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The main goal of this paper is to improve the model by 
IVitense et al.l (l2010h in several important respects: 



1. Although we do not modify the debiasing a lgorithm and 
stay with the same "true" EKB as defined in IVitense et al.l 
(2010), we re-address the question of how the size distri- 
bution in the present-day EKB that we only know down to 
sizes of ~ 10 km should be extrapolated down to the dust 
sizes. Accordingly, we present the new collisional code runs 
that make different assumptions about the amount of objects 
smaller than ~ 10 km in the current EKB. Besides, these new 
runs include a more realistic material composition (a mixture 
of ice and astrosilicate in equal fractions) and an accurate 
handling of the cross-section of dust grains. This is done in 
Section|2] 

II. We estimate the influence of planets (resonant trapping and 
gravitational scattering) (Section|3]l. 

III. We include the possible effect of ice sublimation (Section|4|i. 

IV. We finally make a detailed comparison of the model with 
the spacecraft in-situ measur ements, including the first re - 
sults of New Horizons (Poppe et al.ll2010l: iHan et al ] l2011h . 
as well as with the therm al emission constraints by COBE 
dGreaves & Wv att' '2010', and references therein). This is 
done in Sections |5]and|6] 

Our results are summarized in Section Q and discussed in 
Section HI 
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2. The dust production model 

2.1. The size distribution in the EKB 

We start with general remarks about the size distribution in 
the EKB and its evolution since the early phases of the 
Solar System formation. Since it is not known how plan- 
etesimals in the solar nebula have formed, their primordial 
size distribution is unclear. In standard coagulation scenarios, 
the bottom-up growth of planetesimals could have resulted 
in a broad size distribution (e.g., iKenvon & Bromlevi 120081) . 
with a more or less constant slope across all the sizes up to 
roughly the size of Pluto. Alternatively, local gravitational in- 
stability in turbulent disks would have produced predominantly 
big (~ 100 km) planetesimals (J phansen et al. 2006, 2007; 
ICuzzi et al. 2008; iMorbidelli et alT l2009). implying a knee in 
the size distribution at such sizes which is indi cated by several 
observations ( Bernstein et al. 2004; Fuentes & HolmanI 120081; 
iFraser & KavelaarsI 12009; Fuentes et al. 2009 ). Next, accord- 
ing to the Nice m odel (Gomes et al. 2005; Le vison et al.ll2008t 
iMorbideilil |20 1 Ol) . the primordial Kuiper belt was compact (be- 
tween 15 and 35 AU) and massive (~ 50 Earth masses). With 
these parameters, the EKBOs with sizes up to hundreds of kilo- 
meters would have been collisionally processed by the time of 
the Late Heavy Bombardment (LHB) in s; 800 Myr from the 
birth of the Solar System. Thus, just before the LHB, the size 
distribution consisted of two parts. The objects smaller than hun- 
dreds of km had a size distribution set by their coUisional evo- 
lution in the early massive EKB, whereas the larger objects re- 
tained a primordial distribution set by their formation process. 
The LHB has then resulted in a dynamic al depletion of the EKB, 
which was obviously size-independent dWvatt et al.ll201 ll) . As a 
result, the entire size distribution must have been pushed down, 
retaining its shape. During the LHB, the EKB has reduce d its 
original mass by a factor of ~ 1000 jLevison et al. I l2008h and 
expanded to its present position. Both the reduction of mass and 
the increase of distance to the Sun have drastically prolonged 
the collisional lifetime of the EKBOs of any given size. As a 
result, during the subsequent 3.8 Gyr only the objects smaller 
than about a hundred of meters in radius (more accurate values 
will be obtained later, see Fig.O experienced full collisional re- 
processing. We conclude that the size distribution in the EKB 
after the LHB, and in the present-day EKB, is likely to con- 
sist of three parts. Objects smaller than a hundred meters must 
currently reside in a collisonal equilibrium, those with radii be- 
tween a hundred meters and hundreds of kilometers inherit the 
collisional steady-state of the massive and compact belt of the 
pre-LHB stage, and the largest EKBOs still retain a primordial 
size distribution from their accretion phase. 

2.2. Setup of the collisionai simuiations 

To obtain the dust distributions in the present-day EKB, which is 
the goal of this paper, we use our collisional code A CT {Analysis 
of Co l lisional Evolution ) (Krivov et al. 2000, 2005; Krivov et al.l 
120061 120081; iLohne et aL. 2008c .MuUer et al., ,2010). ACE simu- 
lates evolution of orbiting and colliding soUds, using a mesh of 
sizes s, pericentric distances q, and eccentricities e of objects as 
phase space variables. It includes the effects of stellar gravity, di- 
rect radiation pressure, Poynting-Robertson force, stellar wind, 
and several collisional outcomes (sticking, rebounding, crater- 
ing, and disruption), and collisional damping. 

If we were able to set an initial size and orbital distribution 
of bodies (i.e., the one after the completion of the LHB) in a rea- 



sonable way, we could simply run the code over 3.8 Gyr to see 
which dust distribution it yields. Setting the initial distribution at 
largest EKBOs, i.e. the third of the three parts of the entire size 
distribution described above, is straightforward. Since the distri- 
bution of these objects remains nearly unaltered since the LHB, 
their initial distribution should be nearly the same as the cur- 
rent one. Accordingly, we populate the ACE bins with the debi- 
ased p opulation of known EKBOs, as described in Vite nse et al] 

dloTol) . 

However, we do not know the second part of the distribution, 
at least for objects between a hundred meters and ~ 10 km where 
no or very few EKBOs have been discovered. Given the lack of 
information on these objects, we choose to extrapolate the size 
distribution to smaller objects with a power law dA^ oc ■? ds. 
The slope q is unknown, so we explore the following possibili- 
ties (thin lines in Fig.[TJ: 

1. Ru n "d" ("Dohnany i extrapolation"). We assume the classi- 
cal D ohnanvH (Il969h law with q — 3.5. This e xtrapolation is 
similar to the one used in I Vitense et al.l (1201 Oh . 

2. Run "f" ("flat extrapolation"). We assume q = 3.0 for s < 
10 km. Run "f"can be treated as a rough proxy for a break 
in the size distribution a t a few tens of kilom eters reported 
in the literature: q =a 1.9 (iFras er & Kavelaars 2009). q ~ 2.0 
dFuentes et al.ll2009l) and q ^ 2.5 (Fuentes & H olman 20081) . 
Keeping in mind that the observed TNOs include several 
populations, and that the knowledge of scattered objects is 
particularly poor (Vitense et al. 2010), we made an addi- 
tional run "fcKfl" identical to "f", but without the scattered 
objects. 

3. Run "n" ("no extrapolation"). Here, we refrain from any ex- 
trapolation, assuming that the system was devoid of smaller 
objects initially. This formally corresponds to <7 ^ -oo. 

To complete specification of the initial conditions for ACE 
simulations, we have yet to set the orbital distributions of the ob- 
jects with s < 10 km. On the absense of relevant observational 
data, we simply assume that these objects inherit the pericen- 
tric distance and the eccentricity from their parent bodies. This 
means that, for every object that resides in a bin {si,qj,ek}, the 
bins {si,qj,ek} are populated with (s//i,)'"^ (I < i) objects (as- 
suming logarithmic size bins). 

As a minimum grain radius, we chose 0.4 ^m and set size 
ratios of the adjacent bins of 1.5 for dust sizes and 2.3 for the 
largest TNOs. To cover the heliocentric distances from 4 AU to 
400 AU we used a logarithmically spaced pericenter grid with 21 
bins as well as a linearly spaced eccentricity grid between -1.5 
and 1 .5. Note that negative eccentricities correspond to "anoma- 
lous" hyperbolic orbits, which are open outward from the star 
and are attained by smallest dust grains with a radiation pres- 
sure to gravity ratio > 1 (iKrivov et al.ll2006h . To make sure 
that this, rather coarse, grid yields sufficiently accurate results, 
we made another "n" run with a finer, more extended grid with 
a minimum grain radius of 0.3 \im with size ratios of the adja- 
cent bins of 1.25 for dust sizes and 1.58 for the largest TNOs, 
41 pericenter bins and eccentricity bins between -5 and 5. We 
found that our coarse grid leads to almost the same results as the 
fine grid model. 

As material, we assume d a mixture of 50% ice dWarrenl 
1984) and 50% astrosihcate (^L aor & Drain3ll993h wifli a bulk 
density of 2.35 gcm"^. The optical constants of the mixture were 
computed with the Bruggeman mixing rule and the absorption 
coefficients with a standard Mie algorithm. The values of other 
para meters, e.g. the critic al fragmentation energy, were the same 
as in I Vitense etaP (l2010l) . 
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Fig. 1. Size distributions in different collisional runs. Thin lines 
are initial distributions, while thick ones correspond to an ad- 
vanced state of collisional evolution. Note that the initial size 
distribution of the "n" run coincides with the debiased popula- 
tion of TNOs. A is the cross-section density per size decade at a 
distance of 40 AU. Note that A - const (i.e. horizontal lines) cor- 
respond to a size distribution with q = 3, where different-sized 
objects equally contribute to the cross-section. The gray shaded 
rectangle is a rough approximation of the particle dust flux given 
by New Horizons translated into the cross-section density and 
distances of the EKB. 



2.3. Results of the collisional simulations 

All the extrapolations described above are rather arbitrary, and 
the last one obviously unrealistic. A natural question is then, 
which of the models, and after which timestep, should deliver 
flie distributions that match the actual distributions of the EKB 
material the best. We start with the integration time. Each of the 
runs was let to go as long as needed to reach a collisional equilib- 
rium at smaller sizes, but not too long in order to preserve the ini- 
tial distribution of larger objects. A boundary between "smaller" 
and "larger" sizes was arbitrarily set to s ~ 1 km. We considered 
"collisional equilibrium" to have been reached, once the shape of 
the size distribution stopped changing. Note that, to meet these 
criteria in the "n" run, we had to let the system evolve much 
longer than the age of the Universe. Of course, this "modeling 
time" should not be misinterpreted as physical time of the EKB 
evolution. This was simply the time needed for the population 
of large bodies to generate sufficient amount of smaller debris 
down to dust sizes. 

The results obtained over the integration interval chosen in 
this way are shown in Figs. [T}{3] with thick lines. These three fig- 
ures show the size distribution, the radial profile of the normal 
gemetrical optical depth, and the collisional lifetime of the ob- 
jects, respectively. We note that at an earlier stage of evolution 
the cross-section density and the normal optical depth would be 
lower, and the lifetime of dust grains longer, while a later stage 
of evolution would lead to more dust and therefore to a higher 
cross-section density and optical depth and reduced lifetime of 
the particles. 

But which of the models, "d", "f", of "n" — if any — 
matches the actual dust distribution in the present-day EKB 
the best? The only way to answer this question is to compute, 
for each of the simulations, the observables and compare them 



Fig. 2. Normal optical depth for the same ACE runs and time 
instants as in Fig. [T] 
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Fig. 3. Lifetimes of dust grains and parent bodies for the same 
ACE runs and time instants as in Fig. [T] Particles below the 
3.8 Gyr line are in a collisional equilibrium after the LHB. A 
steep rise in the lifetime at i ~ 300 m corresponds to the 
strength-gravity transition of the critical disruption energy. 



with in-situ spacecraft measurements and thermal emission con- 
straints. Although an in-depth analysis of the data is deferred to 
Sec.|5j we now take a first quick look. The gray shaded rectangle 
in Fig.[T]is a rough approximation of the dust flux data collected 
by New Horizons, translated into the cross-section density and 
extrapolated to the distance of the classical EKB. A compari- 
son with the evolved curves demonstrates that the "d" run is far 
too dusty. It cannot reach an evolutionary stage which would 
be consistent with the measurements (and with the upper limit 
from the non-detection of the thermal emission). Thus we con- 
clude that a straightforward extrapolation from debiased EKBOs 
to dust sizes has to be ruled out. Therefore we have shown here, 
with a compeletely different type of argument, that a break in 
the size distribution has to be present in the EK B, as found from 
the analysis of TNO observations ([F uentes & HolmanI 120081 : 
IFraser & Kavelaarsll2009t iFuentes et"al.i,2009i) . 
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How about the other runs? Both the "f" and the "n" runs 
are consistent with the observational data; we will confirm this 
in Sec.|5]by a more thorough analysis. Thus — unfortunately — 
we cannot constrain the size distribution of EKBOs more tightly. 
Nor can we say which of the dust distributions, the one of the 
"f" run or the "n" run, can be expected in the EKB, although the 
shape of the curves in these runs is different. (The only common 
feature shared by all the curves is an abrupt drop at a: 0.5 \im, 
which is the limit, below which the grains are swiftly removed 
from the system by radiation pressure.) We now come to an anal- 
ysis of these differences. 

The size distribution in the "d" run, which we rejected due 
to violation of observational constraints, is typical of a collision- 
dominated disk. At all sizes, the dust transport is less effi- 
cient than the collisional grinding, and th e cross-s ection den- 
sity peaks just abo ve the blowout limit (iKrivov et al. 2006; 
iThebault & A ugereau 2007). This is also confirmed by the ra- 
dial profile shown in Fig. |2] The outer slope of a; 1 .2 is in a good 
agreement with an approximate analytic solution for a c ollision- 
dominated disk that predicts a slope of a 1.5 (Stru bbe & Chiand 
[2006). Also, there is a clear decrease of the optical depth to- 
ward the star, caused by collisional elimination of the parti- 
cles. Since we can rule out this extrapolation further results for 
the "d" run are presented but not discussed anymore. As men- 
tioned above, the "d" run is essentially the same as the run 
in|yitense et al. (2010) (made for the expected EKB, with the 
Poynting-Robertson effect included), so for a detailed analysis 
of the "d" run we refer to our previous work. As shown here, 
this run fails to describe the actual present-day EKB in the Solar 
System. Neve rtheless, the results and conclusions presented in 
IVitense et al.l (2010) would still be valid for an EKB analog, in 
which all objects down to kilometer in size are in collisional 
equilibrium. 

The size distribution in the "n" run is different. It shows 
a broad maximum at ~ lOOf^im, which indicates that particles 
smaller than that are transported inward from the dense p art of 
the disk before they are lost to collisions (IWvatt et al.l201 ih . The 
inner part of the radial profile in Fig.|2]is nearly constant, and the 
outer one reveals a steeper slope of a 3.0, a s predicted an alyt- 
icall y for a transport-dominated disk (a 2.5. IStrubbe & Chiang! 
l2006i) . Note that the outer profiles are generated by particles in a 
narrow range of sizes around the blowout limit. The coarse size 
grid in our models therefore limits the accuracy with which we 
can reproduce these slopes. 

The "f" run seems to be intermediate. Although the max- 
imum in the size distribution is broader than in the "d" run, 
it still resembles the curves typical of collision-dominated 
disks. However, the profile of the normal optical depth (Fig. |2]l 
stays nearly constant inside th e main belt, which is typical of 
transport-dominated disks (e.g. IWvatil2005h . 

Fi gures [TH3] also present the results of the additional "fcKfi" 
run, from which we excluded scattered objects as dust parent 
bodies. Figure[T]shows, somewhat unexpectedly, that the results 
of "f" and "fcKfi" runs differ from each other: the dust disk in the 
latter turns out to be transport-dominated, similar to the "n" run. 
The question is why. This is not because dropping the scattered 
objects just reduces the amount of material in the EKB, resulting 
in reduced collisional rates. A test simulation, in which we artifi- 
cially augmented the mass of the classical EKB to the total mass 
of the expected EKB, brought qualitatively the same results as 
the "fcKfi" run. Instead, the answer can be found in the method 
of extrapolation. As explained before, we filled the {s, q, e)-bins 
with our debiased population of EKBOs and extrapolated to- 
wards smaller sizes with a power law into the same (q, e)-bins. 



That means that we transfer the high eccentricities of the large 
scattered objects to all smaller ones. Although higher eccentrici- 
ties do not lead to higher collisional rates (see Krivov et al. 200^ 
discussion after their Eq. (17)), they increase the relative veloc- 
ities, making collisions more disruptive. In the "f" run a large 
amount of s < 10 ^m particles is produced, leading to a higher 
number and cross-section density for these particles, which in 
turn leads to a higher collisional rate and a shorter collisional 
lifetime for larger particles (Fig.[3]l. Without the eccentric orbits 
of scattered objects ("fcKfi" run), the relative velocities are mod- 
erate, collisions are less disruptive and fewer small particles are 
produced. Therefore, destruction of larger grains becomes less 
efficient, leading to a prolonged collisional lifetime. 

The above discussion demonstrates that it remains unclear 
whether the EKB dust disk is transport- or collision-dominated. 
It is most likely that it is either transport-dominated or intermedi- 
ate between a collision- and transport-dominated disk. However, 
in all the runs considered, the inner part of the dust disk (in- 
side the classical EKB) has a nearly contant radial profile of the 
optical depth of Tj^ ~ 1 x 10"' (Fig. |2]i. (For comparison, the 
in-plane optical depth for r > 10 AU is Ty = 1 . . .2 x 10"^.) 
This suggests that collisions in the inner part of the disk can be 
neglected. This justifies that in this section we first simulated 
a completely planet- and sublimation-free EKB and will include 
the effects of planetary scattering and subhmation later, in Sec.|3] 
anda 

Figure [3] shows the mean collisional lifetimes averaged over 
all distances for the same ACE runs at the same time instants. 
Note that the collisional lifetime in the main belt is much shorter 
than the average one because the density there is much higher 
and therefore collisions are more frequent. The horizontal line 
represents a lifetime of 3.8 Gyr which is the time elapsed after 
the LHB. All grains below this line are in a collisional equi- 
librium in the present EKB. For all simulations this size is just 
about (10...60)m. The distribution of all objects larger than 
that equilibrium size was set before the LHB and cannot be con- 
strained with our collisional model. 

3. Influence of planets 

Giant planets interact gravitationally with dust in the outer 
Solar System. On the one hand, the grains drifting in ward by 
the Poynting-Robertson (P-R) drag (Burns et al.l ll979l) can be 
captured by planets int o outer mean-motion resonances (e.g . 
Liou & Zook 1999; Mo ro-Martm & Malhotral2002ll2003Ll200l 
iKuchner & St ark 2010). On the other hand, the grains that cross 
the planet's orbit can be scattered. Both effects are able to modify 
the size and spatial distibution of dust in the disk. In this section, 
we investigate the efficiency of capturing and scattering. 

3.1. Resonant trapping 

iMustm & WvattI (1201 ih developed a general formalism to cal- 
culate the capture probability of a particle into the first- and 
second-order resonances with a planet. Their theory is valid 
for any convergent differential migration of the particle and the 
planet (for instance, if the particle is drifting inward and the 
planet is migrating outward). Their results are presented in terms 
of the generalized momentum J and a dimensionless drift rate B 
in their paper). 

The generalized momentum is related to the orbital eccen- 
tricity of the particle reaching the resonance location, e, while 
the dimensionless drift rate B can be expressed through the dif- 
ferential change rate of the particle's semimajor axis, flies, and 
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the semimajor axis itself, Ares- In what follows, we make esti- 
mates for the 3 : 2 resonance with Neptune. Using Eqs. (3) and 
(4) of Must ill & Wvatt ( 2011 ), we find the following conversion 
relations: 



J = 5893.36 



-2/3 



B = -0.818921 pi] 



\OTe / V AU Ores 1 AU Myr 



(1) 



(2) 



where m® and denote the masses of Earth and Neptune, re- 
spectively, and On is the semimajor axis of the Neptune orbit. 

We now assume that Ores is caused by P-R drag 
(IWvatt & Whippld[T950l) 

, BGM^ 2 + 3e^ 



(1 



= -815 



y6 



- e2)3/2 
2 + 3e^ AU 



flres[AU](l-e2)3/2 Myr' 



(3) 



where the pref actor 1.3 accounts for the enhancement of P-R 
drag by solar wind drag dBurns et alj 19791) and /3 is the radiation 
pressure to gravity ratio for the particle. The enhancement by 
solar wind drag is included in all further analysis but we will 
call it just P-R drag for brevity. The y6-ratio not only controls 
the drift rate, it also reduces the effective solar mass felt by the 
particle by a factor of (1 -/?). This affects the resonance location, 
so that Ores rcads: 



Inserting Eqs. (O and (|4]l into Eq. (|2]i, the latter takes the form 

j6 2 + 3e2 



B = 1.6 



(l-/?)2/3(l-e2)3/2- 



(5) 



Using the capture prob abilities as functions of J and B 
from lMustill & Wvat3 (1201 ih and applying Eqs. ^ and we 
computed the probabilities as functions of e and /3 (or equiva- 
lently, particle radius s). The results for the 3:2 resonance with 
Neptune are shown in Fig. |4] Although capturing for grains 
s > 0.6 ^im and e < 0.03 seems unavoidable, it is not obvious, 
what is the fraction of particles of those sizes that will actually 
have such low eccentricities. The reason is that small particles, 
when released from parent bodies in nearly-circular orbits, are 
sent by radiation pressure into large and highly-eccentric orbits. 
Subsequently, drag forces reduce the semimajor axes and eccen- 
tricities of the grains. Yet, it is not clear how low the eccentric- 
ities will be by the time when the grains will have reached the 
resonance location. 

To find this out, we first consider parent bodies with elements 
flp and gp and compute the initial semimajor axis and the ec- 
centricity ei of a grain up on rel ease. To this end, we use Eqs. 
(19)-(20) of lKrivov et al.l (l2006h . in which we neglect the mass 
of the projectile compared to the mass of the target, i.e. the par- 
ent body, and assume that ejection occurs at the pericenter of the 
parent body orbit: 



Gi — a 



(l-/j)(l-gp) 
" l-ep-2yS 

1-/3- 



(6) 



(7) 



Subsequently, the P-R drag will decrease a, and e, . Denoting by 
ef the final eccentricity - i.e. the one the grain will have at the 
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Fig. 4. Capture probability of a single particle with given /3 and 
e at the location of the 3:2 resonance with Neptune. 



location of a reso nance, Qres - cif - and usi ng the dependence of 
fl on e as given in I Wvatt & Whippl3 (Il950h 



de _ 5 e(l - e^) 
da~ 2a 2 + 3e^ 



leads to 



4/5 



1 +, 



Of 



1 - 



(8) 



(9) 



As an example, a plutino with Gp ^ 39 AU and gp =0.1 will 
release a y6 = 0.3 particle into an orbit with a,- = 82 AU and 
e, - 0.57. The 3 : 2 resonance with Neptune for this particle 
is located at Oies ~ 35 AU. At that location, the grain will have 



0.13. 



With Eqs. (6) - (9) and the data of iMustiU & WvattI (1201 ih 
it is possible to calculate the capture probablility for each reso- 
nance and particle size for given Op and gp. As a word of caution, 
we note that the actual dust dynamics can be more complicated. 
One complication is that the initial semimajor axis (Eq. |6l) for 
sufficiently small particles is often so large that the grain has to 
pass several other resonances before it reaches the 3 : 2 one. At 
these resonances, particles with high migration rates and small 
eccentricities will experience an eccentricity jump when they are 
not captured. As a result, our model will underestimate the final 
eccentricity at the 3 : 2 resonance and so overestimate the cap- 
ture probability. Slow migration rates and small eccentricities 
will result in the opposite effect — an eccentricity decrease and 
a probability increase — so an underestimation of the capture 
probabiUty is also possible dMustill & Wvattil20Tlh . A detailed 
modeling of this problem is beyond the scope of this paper. 

Fig-Hshows the probability of capture into the 3:2 resonance 
with Neptune. The probability is the highest for ep = 0, but even 
in that case it does not exceed ^ 20% for dust grains below 2 ^im 
when released from classical EKBOs. Increasing the eccentric- 
ity of the EKBOs and decreasing the grain size reduces the cap- 
turing efficiency. For gp = 0.6, which can be considered repre- 
sentative for scattered disk objects, the capturing probability is 
just a few per cent. For gp = 0.1 (typical classical EKBOs) and 
s ~ I \.im (just above the threshold of the New Horizons dust de- 
tector), the trapping probability is still below 10%. Given these 
results, resonant capturing can be considered unimportant for the 
purposes of this paper and will be neglected. 
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Fig. 5. Capture probability of a /? particle released by a parent body for the 3:2 resonance with Neptune, gp = 0.6 shall represent 
the population of scattered objects. Note that all particles above the horizontal lines have initial eccentricites of e > 1 and will be 
removed from the system (Eq.|7]i. 



3.2. Gravitational scattering 

Since P-R drag continuously decreases the particle's distance 
from the Sun, the grain will eventually reach the orbit of a planet. 
As this happens, the grain can either fall onto the planet, be scat- 
tered, or pass the planet without interaction. In the first two cases 
the particle will be lost. To determine the surviving fraction we 
used a numerical code that calculates the orbital evolution of a 
single particle, taking into account the gravity of one planet and 
the P-R effect. For each y6-value listed in Table [T] (these are the 
same values as used in our collisional simulations) we started 
10.000 particles, with an EKB-l ike a, e, i distr i bution taken from 
the upper panel of Fig. 5 from I Vitense et alj (|2010|) . A particle 
was counted as a survivor as soon as its apocentric distance be- 
came smaller than the pericentric distance of the planet. 

The results are listed in Table [1] for Neptune, Uranus and 
Saturn. As expected, the surviving rate decreases for larger 
grains with lower migration rates. For Neptune and Uranus the 
ejection rate is negligible and will not alter the dust flux signifi- 
cantly (Sec. Is]). However, Saturn ejects nearly a half of the dust 



Table 1. p values and corresponding sizes, masses and surviving 
rates for particles passing Neptune, Uranus and Saturn. Particles 
between 10"^^ g < m < 10"^ g can be m easured by the New 
Horizons dust counter (iHoranvi et alJl2008h . 



p 


s [\im] 


m [g] 




^surv 




^surv 


0.404 


0.65 


2.7 X 10- 


12 


96.8% 


97.4% 


79.7% 


0.259 


0.99 


9.5 X 10- 


12 


93.8% 


95.1% 


66.0% 


0.164 


1.5 


3.4 X 10- 


11 


88.7% 


90.0% 


57.0% 


0.106 


2.3 


1.2 X 10- 


10 


82.2% 


80.1% 


50.3% 


0.070 


3.5 


4.3 X 10- 


10 


78.3% 


73.9% 


47.1% 



grains. As we will see in Sec.|5j Saturn's influence is important 
to explain the in-situ measurements, but all three planets have 
little effect on the thermal emission of the EKB dust (Sec.|6]l. 

As shown in the previous section, the EKB dust disk is 
transport-dominated for small particles, which means that col- 
lisions play a minor role. Therefore, gravitational scattering can 
simply be implemented by multiplying the distribution obtained 
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in the collisional simulation by the surviving rates for the corre- 
sponding particle sizes and distances. 

4. Influence of sublimation 

When drifting inward, dust grains will not only suffer interac- 
tion with planets, but they will also be heated up due to the de- 
creasing distance to the Sun. Our dust particles are composed of 
"dirty ice" (50% ice and 50 % astrosilicate in volume). Their ic y 
part sublimates at ^ lOOK (!KoMiashi^l.''2008','2009^ 
Since the EKB dust disk is radially optically thin, the temper- 
ature of a dust grain is determined by the energy balance be- 
tween the absorption of incident solar radiation and the thermal 
emission of the grain. We neglect the latent heat of sublimation 
because its contribution is minor (Kobav ashi etanr2008,) . The 
sublimation distance rsubi, where the temperature of a particle 
reaches 100 K, depends on its size. If the particles are larger 
than A/ (In), where A is the peak wavelength of emission, the 
absorption and emission cross-sections are approximately the 
same as the geometrical one, and these particles can be assumed 
to be blackbody radiators. Since for T = 100 K the maximum 
is at /I ~ 30 ^m, this is the case for grains with 5 > 5 ^m. 
Temperatures of smaller particles are obtained by solving the 
thermal balance equation (see, e.g., Krivov et al. 2008). Fig. |6] 
shows the resulting temperatures for different sizes and dis- 
tances, with three isotherms overplotted. The leftmost one cor- 
responds to lOOK. Empirically we can approximate the depen- 
dence of the sublimation distance (in AU) on the size (in mi- 
crometers) by: 



'"subl 



-10.2 sin(0.26.?) -i- 16.85 s < 5.0 \im 
8.0 s > 5.0 |xm 



(10) 



Note that this function does not have a physical meaning and is 
only needed to implement sublimation into our model. 

The outcome of sublimation depends on the structure of icy 
grains. If a single icy particle is an aggregate of small grains, 
each having /3 ^ 0.5, the resulting grains will be blown out and 
therefore no grains should be present inside r^^bi- However, if 
the constituent monomers have /3 ^ 0.5, the number density of 
grains inside rsubi will increase. Since both is inconsistent with 
the dust flux measured by spacecraft, a single icy grain is likely 
to contain a single core of refrac tory material covered with an 
ice mantle (iKobavashi et al.ll2010l) . For our dirty-ice grains sub- 
limation will result in a 100% silicate particle which has a half 
of the volume of the original particle. The radius of the resulting 
particle is simply 

^silicate = VoiSiicy, (11) 



and the mass is given by 



^^silicate 



r, ^Psilicate 

: 0.5 m. 



Picy 



icy 5 



(12) 



with Picy = 2.35 gem ^ being the bulk density of the dirty ice and 
Psiiicate = 3.35 gcm""' of the astrosilicate. The typical sizes and /3- 
values of the particles before and after sublimation, together with 
their sublimation distances, are given in Table |2] 

We now discuss how sublimation affects the distribution of 
dust. The particles born through collisions in the Edgeworth- 
Kuiper belt have eccentricities roughly comparable to their/? val- 
ues (Eq.|7]i. Although damped by P-R drag, their eccentricities in 
the sublimation zone are typically larger than 0.05. Particles with 
e > 0.05 will experience a rapid sublimation without pile-up 



T, [Kl 



50 




Fig. 6. Temperatures of the dirty-ice particles for different dis- 
tances. The solid lines correspond to 100 K, 75 K and 50 K. 
Sublimation occurs at Tsubi = lOOK. Sublimation distance in- 
creases with decreasing size because the emission efficiency of 
small grains is lower, which makes them hotter. 

Table 2. Sizes and yS-values before and after sublimation and 
the corresponding sublimation distances; for particles larger than 
5.0 \im blackbody temperatures are assumed. 



■Sicy [kim] 


/Sicy 


'■subi [AU] 


■^silicate [^im] 


^silicate 


0.425 


0.576 


15.7 


0.337 


0.652 


0.648 


0.404 


15.1 


0.514 


0.428 


0.989 


0.259 


14.3 


0.785 


0.280 


1.51 


0.164 


13.0 


1.20 


0.184 


2.30 


0.070 


11.1 


1.83 


0.121 


3.51 


0.046 


8.8 


2.79 


0.079 


>5.0 




8.0 







and d ust ring formation (IKobavashi et al.ll2009l) : (cf. lBurns et al.l 
1979, their Fig. 8). Next, although sublimation in our model does 
not eliminate the particles and thus preserves their number, it 
reduces their spatial number density. This is because the num- 
ber density of particles is inversely proportional to their drift 
rates in the steady state. Because a cc p, the increase of /3 due 
to sublimation lessens the number density of particles. Based 
on Table |2] the change is estimated to be only about 20%, see 
Fig. Q below. However, with in-situ dust detectors measuring 
only grains above a certain threshold, the observable dust flux 
decreases more strongly. 

If we assume that the orbital changes due to the change in 
size and therefore changing interaction with the stellar radiation 
are small, we can implement sublimation into our collisional re- 
sults the same way as gravitational scattering by simply correct- 
ing sizes and cross-section- and mass density for the affected 
bins. Since planetary scattering and sublimation are independent 
processes, the order of implementation does not matter. 
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5. Comparison with spacecraft measurements 

The Student Dust Counter on-board the New Horizons space- 
craft is capable of detecting impacts of grains with 10 '^g < 
m < 10"^ g and can distinguish grain rnasses apart by a factor 
of 2 between 0.5 \im < s < 5 [im (fHoranvi et alj 12008 ). The 
first results from iPoppe et alJ (1201 Ol) i ndicate pa r ticle fl uxes up 
to 1.56 X 10-4 m-^s"'. The results of iHan et all (1201 ll) show a 
slight increase of the flux for r > 15 AU. The particle flux can 
be calculated via 



mnVrei d(ln m) 



(13) 



with m being the mass of the particle, n the number density 
per logarithmic mass and Viei the relative velocity between the 
spacecraft and the particle. The first two values are a direct out- 
put of our simulation. The relative velocity was assumed to be 
v'rei = 15.54 kms"', according to the official New Horizons web 
pag43- Based on the results of Secs.|2]3 we calculated the dust 
fluxes for the EKB dust disk unaffected by planets and sublima- 
tion, the one with planets and the one with planets and sublima- 
tion. Although the separate contributions of planets and sublima- 
tion are rather low, their combination can alter the dust flux up to 
a factor of three, whereby Saturn plays the most important role. 
In Fig.|7]the results of run "f" are shown. The right evolutionary 
state of the simulation (i.e., timestep) was chosen in the follow- 
ing way. As seen in Fig.|7] the black solid line can be assumed 
to be a constant for r < 20 A U. Taking thi s assum ptio n we fitted 
the N ew Horizons data from iPoppe et"al] (1201 Oh and iHan et all 



201 ll) by a constant line to F^^^^^ « 3 X lO"'* m'^s"' and searched 
for the timestep w hich agrees wi t h the m odel the best. 

According to iGurnett et al.l (Il997t) . the Voyager 1 and 2 
plasma wave instruments, that acted as "chance" dust detectors, 
have a mass threshold of m > 1.2 x 10 g, which is one order 
of magnitude higher than for the New Horizons dust counter 
Accordingly, we rescaled the Voyager data to the New Horizons 
threshold, with a power law slope of <7 = - 1 obtained in our sim- 
ulation for the corresponding masses (Fig.lTJ at 40 AU. Since the 
instruments aboard Voyager I and II were neither designed to de- 
tect dust impacts nor calibrated for this purpose and traversed the 
outer Solar System in highly inclined orbits, their dust measure- 
ments should be compared with our model with great caution. 

Simulations f, fcKB, and n were treated the same way. Since 
for the particle sizes in question {s < 5 (xm) all modeled disks are 
transport-dominated, the results do not differ much from each 
other and lead approximately to the same fits as for the "f" run. 
Therefore these results are not shown in Fig.|2] 

6. Thermal emission constraints 

To calculate thermal emission of dust in the EKB, we computed 
the photospheric spectrum of the Sun using the NextGen mod- 
els (iHauschildt et al.l [19991) . The eq uilibrium dus t temp eratures 
were obtained by the procedure of iKrivov et a n (I2008h . As in 
the coUisional simulations, we adopted the "dirty ice" consisting 
of equal volume fracti ons of ice (Warren 1984) and astrosilicate 
dLaor & Drainel[T993h . As explained in Sec. |4]the sublimation 
distance depends on the particle size. Therefore we divided the 
EKB into 6 sub-rings to handle the different emission properties 
of the dirty ice and pure astrosilicate (TableO. 

To place our EKB in the context of extrasolar debris disks, 
we now consider the EKB dust disk, as if it were viewed from 



' http://pIuto.jhuapl.edu/mission/whereis.nh.php (Last accessed on 
2 September 2011) 
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Fig. 7. Simulated particle flux in comparison with in-situ mea- 
surements by New Horizons and Voyager 1 and 2. The solid 
black line takes into account planets and s ublimation and rep - 
resents ou r best fi t to t he data taken from Poppe et alJ (l2010l) . 
iHan et alJ (1201 ih and iGurnett et all (Il997 i). The dotted and 
dashed lines show the dust flux for the unperturbed EKB and 
after planetary scattering without sublimation, respectively. 

Table 3. The EKB divided into 6 sub-rings: material composi- 
tion together with the sizes and distances for which that compo- 
sition was adopted. 



ring # 


material 


s [\.im] 






distance [AU] 


1 


astrosilicate 


0.425 < 


s 


< oo 


< r < 8 


2 


dirty ice 


3.51 < 


s 


< oo 


8 < r < 12 


3 


astrosilicate 


0.425 < 


s 


< 3.51 


8 < r < 12 


4 


dirty ice 


1.51 < 


s 


< oo 


12 < r < 16 


5 


astrosilicate 


0.425 < 


s 


< 1.51 


12 < r < 16 


6 


dirty ice 


0.425 < 


s 


< oo 


16 < r < oo 



10" 



10" 




run f 

run f w planets 

run f w planets & subl 



10 




1000 



Fig. 8. Spectral energy distribution of the EKB including planets 
and sublimation (solid black line), without sublimation (dashed 
red line), and for an unperturbed EKB (dottted blue line). All the 
curves are based on the same run ("f") and the same time instant 
as in Fig.|2l 
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outside. The final spectral energy distribution (SED) of the EKB 
dust disk, as seen from a reference distance of 10 pc, is shown 
in Fig. [8] The influence of planets and sublimation does not al- 
ter the shape, peak position, and height significantly. The SED, 
coiTected for planets and sublimation, peaks at 40 - 50 [im with 
a maximum thermal emission flux of ^ 0.5 mJy, which amounts 
to * 0.5% of the photospheric flux at that wavelength. The pre- 
dicted flux drops to a: 0.4 mJy at 10 [im and to a: 0.2 mJy at 
100 [im. The fractional luminosity of the modeled EKB dust 
disk, after applying corrections for planets and sublimation, is 
fd = 1.2 X 10-^. 

Our results are consistent with the upper limit, placed by 
non-detection of the EKB dust emission at 70 ^im with the COBE 
spacecraft. T hat limit amounts to 1 + 0.5% of the solar photo- 
spheric flux (Gre aves & Wvattll2010l) and is shown by the dark 
gray box in Fig.|8] 

Would EKB analogs around nearby stars be detectable, for 
example, with the PACS instr ument CPoglitsch et al..2010 ) of the 
Herschel Space Observatory (iPilbratt et al.ll2O10l) ? The sensitiv- 
ity of the PACS instrument at 70 \im is 4.7 mJy in 1 hour inte- 
gration time at a 5cr uncertainty leveS This is about a factor of 
10 above the calculated SED flux of 0.4 mJy at 70 i^im. This fac- 
tor would further increase when taking into account additional 
background noise and photospheric flux uncertainties. We con- 
clude that the detection of an exact analog of the EKB with the 
presen t-day instruments is impossible. An apparent contradic- 
tion to lVitense et al.l (1201 Ol) . who concluded that Herschel/PACS 
should be able to detect an ^ 2Mekb analog, traces back mainly 
to a different extrapolation method from parent bodies to small- 
est grains. 



7. Conclusions 

The purpose of this paper was to develop a self-consistent model 
of the EKB debris disk. To accomplish this task we used the 
debias ed population of EKBOs as described in ' Vitense et alj 
(I2OIOI) . Treating this population as dust parent bodies, we gen- 
erated their dust disk with our coUisional code. We draw the fol- 
lowing conclusions: 

1 . We have shown that sub-kilometer-sized EKBOs largely de- 
termine the amount and distribution of dust in the outer Solar 
System. However, these are far too small to be directly de- 
tected at present in TNO surveys, and their properties cannot 
be accessed by coUisional modeling, because they are not in 
a coUisional equilibrium. Therefore, an extrapolation from 
observable TNOs towards smaller sizes is necessary. 

2. A straightforward extrapolation for the yet unknown objects 
(s < 10 km) with a classical Dohnanyi law can be ruled out. 
In that case, the amount of dust would be so large that its 
thermal emission would have been detected by the COBE 
spacecraft. Therefore, the distribution of these objects should 
be flatter. In other words, a break in the size distribution at 
several tens of kilometers has to be present. 

3. Different extrapolation methods which are consistent with 
the measurements reveal the EKB either as a transport- 
dominated debris disk or to be intermediate between 
the collision-dominated and transport-dominated regimes. 
Depending on the extrapolation method, we found the 
present-day EKB to be in coUisional equilibrium for objects 
s < (10...60)m. 



4. Using the results of iMustiU & WvattI (1201 ih . we estimated 
the effect of resonance trapping of planets. The capturing 
rate of the dust grains that are either detectable with in- 
situ measurements by spacecraft or contribute to measurable 
thermal emission turned out to be < 10% in most cases and 
not to exceed < 20% even for the largest grains considered. 
Accordingly, resonance trapping should have a negligible ef- 
fect on dust impact rates and dust thermal emission, given the 
typical accuracy of the dust measurements. 

5. Gravitational scattering of dust grains by planets was inves- 
tigated numerically. Scattering can modify the particle flux 
in the Saturn-Uranus region (8 AU < r < 15 AU) by about 
a factor of two and has little effect on thermal emission of 
dust. 

6. Likewise, sublimation can reduce the particle flux by approx- 
imately a factor of two and does not affect the thermal emis- 
sion fluxes perceptibly. 

7. We calibrated our model with the in-situ measurements of 
the N ew Horizons dust counter (iPoppe et alJ2010l : [Han et alJ 
l20Tlh by fitting our results to the data points and can re- 
produce the nearly constant particle flux of 3 x 10""* m^^s"' . 
The corresponding production rate of dust inside the EKB 
amounts t o 2 x 10^gs~', consisten t with previous esti- 
mates (e.g.lYamamoto & Mukai|[T99l iLandgraf et alj|2002l: 



iHan et alJl201 ij). In a steady-state coUisional cascade (which 
we assume), the "dust production rate" is the same as the 
"dust loss rate". Thus the result means that 2 tons of dust 
per second leave the system by inward transport and through 
ejection as blowout grains. 
8. The spectral energy distribution of an EKB analog, seen 
from a distance of lOpc, would peak at 40 - 50 ^im with 
a maximum flux of 0.5 mJy. This is consistent with the up- 
per limit, placed by non-detection of thermal emission from 
the EKB dust as it would be viewed from outside at 70 (im 
by the COBE spacecraft. The fractional luminosity of the 
EKB was calculated to be fj = 1.2 x 10"^. The in-plane op- 
tical depth for r > 10 AU is set by our model to 2 x 10"^. 
Although the Herschel/PACS instrument successfully de- 
tects debris disks at similar fr actional luminosity levels as the 
EKB (Eiroa et al. 2010, 201 1), all these are larger and there- 
fore colder. Their thermal emission peaks at wavelengths 
longward of 100 (xm, where the stellar photosphere is dim- 
mer. The detection of an exact EKB analog even with PACS 
would not be possible. 

8. Discussion 

Like every model, ours rests on many assumptions and is not 
free of uncertainties. Here we discuss some issues. 



^ http: //herschel. esac. esa. int/Docs/PACS/html/ch03s05. html#sec- 
photo- sensitivity (Last accessed on 24 November 20 11) 



1. Material composition. Compared to IVitense et alj ( 1201 Ol) . 
who applied geometric optics in calculating the radiation 
pressure, we now consistently used a more realistic mate- 
rial composition in both coUisional and thermal emission 
calculations. Nevertheless we assumed many parameters of 
solids — bulk density, shape, porosity, tensile strength, and 
others — to be the same across the entire size range, from 
Pluto-sized TNOs down to dust. This assumption is obvi- 
ously unrealistic. 

2. The role of sub-kilometer-sized EKBOs. Since little is 
known about EKBOs smaller than a few tens of kilometers, 
but these largely control the amount and distribution of dust 
in the outer Solar System, an extrapolation from observable 
TNOs towards smaller sizes is necessary. The question is 
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what kind of extrapolation is reasonable. If the parent bod- 
ies inherit their orbital elements to their children and grand- 
children, then the EKB should comprise a huge amount of 
meter- and sub-kilometer-sized objects in highly eccentric 
orbits, stemming from scattered EKBOs. This would make 
collisions more disruptive and alter the size distribution of 
dust. The resulting size distribution would be dominated by 
the smallest dust grains, just above the radiation pressure 
blowout limit. If, in contrast, the meter- and sub-kilometer- 
sized objects have moderate eccentricities, the peak of the 
cross-section in the size distribution would be broader and 
shifted to larger grains. To distinguish between these pos- 
sibilities, one needs more information about the amount 
and distribution of sub-kilometer objects. Accurate measure- 
ments of sizes and orbital elements of dust grains in the outer 
Solar System in the future would also help. 

3. Break in the size distribution of EKBOs. Our model shows 
that a break in the size distribution at tens of kilimeters, as 
reported is the recent literature, is necessary. Otherwise the 
EKB dust disk would be too dusty, violating the available ob- 
servational constraints. If such a break is present in other de- 
bris disks as well, then the total mass of parent bodies should 
be higher that that usually inferred in the debris disks studies. 

4. Planetary scattering and sublimation. If planetary scatter- 
ing and/or sublimation is more efficient than assumed in 
our model, the amount of dust grains that reach the Saturn- 
Uranus region of the Solar System would be smaller. To stay 
consistent with the in-situ measurements, one would have 
to compensate higher scattering rates and/or more efficient 
sublimation by higher dust production rates in the EKB. 
However, this would lead to an higher thermal emission flux, 
which would contradict the non-detection of thermal emis- 
sion by COBE. 

More and deeper TNO surveys, including good measure- 
ments of their orbital elements (in particular the eccentricity), 
would help to improve the extrapolation method and therefore 
our model, resulting in tighter constraints on the dust distribution 
and more accurate predictions for the upcoming dust flux mea- 
surements. Better constraints on the population of sub-kilometer 
objects, which could be expecte d, for instance, from the stellar 
occultation m ethod (e.g.. lLiu et a l. 2008; Schlichting et al. 2009; 
iBianco etal ][201fl), would also be of great help. Of course, the 
most promising way to improve the model of the EKB dust is to 
use direct observations of dust. In particular, a size distribution 
of the impacted grains on the New Horizons dust counter would 
be very helpful as well as new thermal emi ssion contraints that 
could be expected from the Planck mission (lAde et al.ll20TTh . 
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